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Abstract. In this paper, we study the initial-value problem for two first order 
systems in non-conservative form. The first system arises in elastodynamics 
and belongs to the class of strictly hyperbolic, genuinely nonlinear systems. 
The second system has repeated eigenvalues and an incomplete set of right 
eigenvectors. Solutions to such systems are expected to develop singular con- 
centrations. Existence of singular solutions to both the systems have been 
shown using the method of weak asymptotics. The second system has been 
shown to develop singular concentrations even from Riemann-type initial data. 
The first system differing from the second in having an extra term containing 
a positive constant k, the solution constructed for the first system have been 
shown to converge to the solution of the second as k tends to 0. 



1. Introduction 

The initial- value problem for the first-order quasilinear hyperbolic system 

du du da 

dt dx dx /, ,| 

da da 2 du 

"TTT + U- k — = 

dt dx dx 

(in the domain f2 = {(x,t) : — oo < x < oo,t > 0}) arising in applications in elas- 
todynamics, has been well studied (see [2], [5], [6], [8]). Here A: is a positive constant. 
It is a strictly hyperbolic system having two real distinct eigenvalues given by 

Ai (it, a) = u — k, A2 (u, a) = u + k 

with the corresponding right eigenvectors 

E l {u,a)= ( * V E 2 (u,a)= ( \ 

Now letting k — > 0, we see that the eigenvalues X\(u, a) and A^U, a) tend to 
coincide. In particular, taking k = in (jl.ll) we arrive at the system 

du du da 

dt U dx dx ~ ' (1 2) 

da da 

m +u dx- = Q - 
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which has repeated eigenvalues Xi(u, cr) = A2 (u,a) = u and an incomplete set of 



In [llj . a class of 2x2 systems in conservative form having an incomplete set of 
eigenvectors everywhere has been considered. These systems exhibit development 
of singular concentrations. We expect a similar kind of development of singular 
concentration for the system (|1.2[) . But the analysis in pj] cannot be applied di- 
rectly due to the following two reasons: 

1. The system (|1.2[) is in nonconservative form and hence we need to give a suitable 
meaning to the nonconservative products which in general lead to different solu- 
tions depending upon the meaning attached (see [3].[9].^0]1. 

2. One of the assumptions that has been used in [11] is that the eigenvalue should 
have vanishing directional derivative along a right eigenvector. But the eigenvalue 
u in this case has a nonvanishing directional derivative along any right eigenvector. 
Thus the assumption is not satisfied. 

Another reason to expect singular solutions for the system (|1.2|) comes from 
studying the behaviour as k — > of the shock and rarefaction curves obtained for 
the Riemann problem for the system (|1.1[) . In [5J, the Riemann problem for the 
system (jl.ip has been studied using Volpert's product. Starting with Riemann type 
initial data 



the shock curves Si(ul,(Tl), S^Ul, o~l) an d the rarefaction curves R\(ul, (Jl), 
^(uljPl) can be written down in the u — a plane as in [6j: 

Ri(ul,ctl) ■ cr = &L + k(u - ul), u>u Ll 

R2(ul,ctl)- a = a L - k(u - u l ), u>u l , 

Si{u L ,a L ) : a = a L + k(u - u L ), u<u L , 

S2{ul,°l) ■ cr = cr L - k(u- u L ), u<u L . 

Now as k — > 0, it seems from the above curves that we cannot have a jump in a 
and therefore the Riemann problem might not be solvable using Volpert's product. 
That this is indeed true when k = 0, that is, for the system (|1.2[) follows from the 



For the solution of the Riemann problem using Volpert's product the following 
relations have to be satisfied (see [3]) 



where s denotes the speed of the discontinuity and [w] denotes the jump in the 
function w across the discontinuity. It follows from the above relations that [cr] = 
thus proving what we expected. 

This observation again makes us suspect that there might be development of sin- 
gular concentrations even if we start with Riemann type initial data. 

In what follows our main aim is to construct generalised 8— shock wave solutions for 





following: 



«N + [y] - M = 0, - s[a] + -iu L + u R )[a] = 
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the systems and (|1.2p and examine the role played by the constant k. We use 
the method of weak asymptotics (see [l],[I],[7]) as the main tool. Depending upon 
the cases, we shall sometimes use the complex-valued weak asymptotic solutions 
(see [7]) for our construction. 

The plan of the paper is as follows: In Section 2, the notions of the weak asymptotic 
solutions, generalised S— shock wave type solutions are recalled and a brief sketch 
of the method of weak asymptotics is given. In Section 3, a few weak asymp- 
totic expansions are proved which are crucially used in the construction of singular 
solutions. In Section 4, existence of generalised 8— shock wave type solutions for 
the systems (jl.ljl and (|1.2p is proved. The role played by k is examined and the 
development of singular solutions from Riemann type initial data is exhibited. 



2. A BRIEF DISCUSSION ON THE METHOD OF WEAK ASYMPTOTICS 

In this section, we recall the definitions of weak asymptotic solutions and gen- 
eralised S— shock wave type solutions and discuss the method of weak asymptotics 
(see[5,@,lZI)- 

Let T> and T>' denote the space of smooth functions of compact support and the 
space of distributions respectively. 

Let Oj)'(e a ) denote the collection of distributions f(x,t,e) € T>'(M) such that for 
any test function <p{x) € T>(M) the estimate 

holds and is uniform with respect to t. The relation ox>'{e a ) is interpreted similarly. 

Definition 2.1. ([7]) A pair of smooth complex- valued (real- valued) functions 
(u(x, t, e), <r(x, t, e)) is called a weak asymptotic solution of the system (jl.ip with 
the initial data (u(x, 0), o~(x, 0)) if 

du(x,t,e) du(x,t,e) da(x,t,e) 
___+„(*, t , e) __ __ =0w{ll 

da(x,t,e) da(x,t,e) 2 du(x,t,e) 

____ + <Xt t, e)—^— - k = o v ,(l), (2.1) 

u(x, 0,e) - u(x,0) = ox)'(l), 
a(x,0,e) - cr(x,0) = ov(l), e 0. 

Similarly, for the system (11.21) we define 

Definition 2. 2. ([7]) A pair of smooth complex-valued (real-valued) functions 
(u(x, t, e), o~(x, t, e)) is called a weak asymptotic solution of the system (|1.2p with 
the initial data (u(x, 0), cr(x, 0)) if 

du(x,t,e) du(x,t,e) da(x,t,e) 
___ + „(«, t , e) __ __ =0p ,(i), 

da(x, t, e) , .da{x,t,e) 

g t + *> £ ) fe = oc(l). (2.2) 

u(a;, 0, e) — 0) = ox)'(l), 

cr(x, 0, e) - a(x, 0) = ov> (1), e -> 0. 
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Definition 2.3. A pair of real- valued distributions (u(x,t),a(x,t)) £ C(R+; V(M.)) 
is called a generalised solution of the systems or (II. 2[) if it is the weak limit 
(limit in the sense of distributions) of a weak asymptotic solution (u{x, f , e), a{x, f , e)) 
as e — > 0. 



2.1. The method of weak asymptotics. Instead of writing it down schemati- 
cally (see [1]), here we shall try to illustrate the method of weak asymptotics by 
applying it to the systems (jl.ljl or (|1.2[) with initial data of the form 

u(x, 0) = u + uiH(-x), 

(2 3) 

<t{x,0) = a + o-iH(-x) + e°5{x). V ' J 

where uo,ui,<7o,<7i, e° are constants. In this case we seek solutions in the form of 
the singular ansatz 

u(x,t) = u + UiH(-x + 0(f)), 

a(x,t) = a + aiH{-x + 0(f)) + e(t)S(x - <f>(t)) 

Such a solution will be called a generalised 5— shock wave type solution for the 
systems (JOJ or (fl~2|) . 

The first step is to regularize the Hcavisidc functions and delta distributions and add 
correction terms so as to form a smooth ansatz for the weak asymptotic solutions 
in the form: 

u(x, t, e) = u Q + uiH u (-x + <f>{t),e) + R u {x, t, e), ^ 

a(x, t, e) = (T + a\H. a [—x + <f>(t), e) + e(t)J(x — 0(f), e) + R a {x, t, e). 

Here H Ul H a are regularizations of the Heavisidc function, 6(x, e) is a regularization 
of the delta function and R Ul R a are correction terms chosen so as to satisfy the 
conditions: 

Ri{x,t,e) = ox>'(l), ^ = op'(l), e ->• 0, i = u,a. 

The next step is to prove the existence of a weak asymptotic solution by substituting 
the smooth ansatz in place of u and a in the left hand side of the systems and finding 
out the functions <j){t) and e(t). 

Once this step is completed, we obtain a generalised 6— shock wave type solution of 
the form (|2.4p as the distributional limit of the weak asymptotic solution as e — > 0. 

3. The regularizations and weak asymptotic expansions 

In this section, we describe the regularizations and correction term to be used in 
the next section and prove a few weak asymptotic expansions for various products 
containing them. 

As mentioned in Section 2, we seek generalised 5 —shock wave type solutions of the 
form 

u(x, t) = Uq + U\H{—X + 0(f)), 

a(x,t) = o-q + cri H(-x + 0(f)) + e(t)S(x - 0(f)) J ' 

for the systems (jl.l[) and (|1.2[) with initial data of the form 
u(x, 0) = Uq + uxH(-x), 

(3-2) 

a(x, 0) = (t + aiH (-x) + e°S(x). 
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Here <f>(t) and e(t) are smooth functions to be determined. To apply the method 
of weak asymptotics, we suggest a smooth ansatz of the form (we do not take any 
correction term in the expression for o~(x, t, e)) 

u(x, t, e) = u + uiH u {—x + <j>(t),e) +p(t)R(x — (p{t), e), 
a(x, t, e) = cr + <J\H a {—x + <j>{t),e) + e(t)S(x — (f>(t),e). 

Here H u (.,e), H a (.,e) are regularizations of the Heaviside function, S(.,e) is the 
regularization of the delta distribution and p{i)R(.,e) is a correction term where 
p(t) is a smooth function (complex-valued or real- valued) to be chosen afterwards. 
The regularizations of the delta distribution and the correction term used in the 
proofs are same as in [7] . The regularization of the Heaviside function has a subtle 
change though. 

Let to : R — ¥ R be a non-negative, smooth, even function with support in (—1, 1) 
and satisfying 



uj{x)dx 



Let wq = j^^ 2 ( x )dx and let 



1 ,x — 2e N „, , 1 ,x + 2e, 

R(x,t,e) = -=w{ ), 6(x,e) = ~w( ). (3.4) 

ye e e e 

Since by definition, R is independent of t, we henceforth denote it by R(x, e). 
An easy calculation using test functions then shows that R(x, e) = od> (e) and since 
R is independent of t, we also have dR g^ — ox><(e). Thus it satisfies the criteria 
for being a correction term as stated in Section 2. 
Also it can be easily checked that 

™&3- =0iy (e), R 2 (x,e)=oj Q S(x)+ov>(e), R(x,e)^^- = )y Q 8' (x) + op, (e) 
and that 

8(x, e) = 6{x) + ow(e), = S'(x) + o v \e). 

ox 

An important observation we need to make at this point is that the supports of 
R(x, e) and 5(x, e) are disjoint. Therefore we have 

dS(x e) 

R(x, e)8(x, e) = 0, R(x, e) — = 0. 

ox 

Letc=(|-^). 

Let us define the smooth function H(x, e) as follows: 




H(x,e) = {c, -3e<x<3e (3.5) 



and is continued smoothly in the regions (— 4e, — 3e) and (3e,4e). 

We take H u (x, e) = H a (x, e) = H(x, e). Again a little bit of calculation shows that 

H(x,e) = H(x)+o v ,(e), f) = 5(x)+o jy {e), H(x,e) ^^ - ls(x)+o v ,(e) 
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Since the supports of i?(x, e) and 5(x, e) are contained in (— 3e, 3e), it again follows 
that 

dR{x,e) dR{x,e) 
HM ^x- = C ^x- =OV ' { ^ 
BH(x e) 

R(x,e) ±£ '- = 0.R(x,t,e)=o v ,(e), 

x dS(x,e) , 
H(x,e) — — — = cb (x) +o V '(e). 

From the above discussions we then have the following lemma 

Lemma 3.1. Choosing the regularizations and corrections as in (|3.4I) and (13.51) we 
have the following weak asymptotic expansions: 

R(x, e) = o v > (1), — — — = o v > (1), 

R 2 (x, e) = uj S(x) + ox)'(l), 

dR^e± = 1 ^ + oc 
ox 2 

S(x, e) = S(x) + o v > (1), = S'(x) + OT> , (1), 

R(x, e)6(x, e) = 0, R(x, e)^|^ = 0, (3.6) 

H(x, e) = H(x) + o v ,(l), dH( f' e) = S(x) + o v ,(l), 

ox 

TT , , 8H(x, e) 1 . ,„ . 

H(x, e) ^ = -5{x) + o v ,{l), 

8R(x,e) dH(x,e) 
H(x,e) — — — =023/(1), R(x,e) — — — = o v ,{l), 

H(x,e)^^- = cS'(x) + o v ,(l), e^O. 
ox 



4. Generalised Delta-shock wave type solutions via construction of 

weak asymptotic solutions 

In this section, we construct 5— shock wave type solutions for the systems (|1.1[) 
and (jl.2l) using the method of weak asymptotics. 

For the rest of the discussion, we use the convention [v] = vl — vr, where vl,vr 
respectively denote the left and right states of v across the discontinuity. 

Theorem 4.1. For t G [0,oo), the Cauchy problem (jl. II) . (13.21) has a weak asymp- 
totic solution (|3.3I) with cj)(t),e(t) and p(t) given by the relations 

0(t) = ^dzM, e( t ) = fi-jfcv, 

N M i (4.1) 

^p 2 (t)uj - e{t) = 0, 

where ujq * s o. positive constant (defined in Section 3 ). 
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Proof. To begin with, let's recall the form of the smooth ansatz: 

u(x,t, e) = u + uiH u (-x + 4>{t), e) + p(t)R(x - <f>(t),e), 
a(x, t, e) = (T + o\H a {—x + </>(£), e) + e(t)S(x — <p(t),e). 

Then the partial derivatives of u(x, t, e) and o~{x, t, e) are given by 

du(x : t,e) ■ dH u (-x + (j)(t),e) x dR{x - (j>(t) , e) 
— = u x cj){t) ^ \-p(t)R{x - (p{t),e) -p{t)(/>[t) — 

du(x,t,e) dH u (—x + 4>(i),c) . . dR(x — <t>(t), e) 
~ = ~ Ul d£ +m cH ' 

= ai ^ t) d Jhk^±MA + ms{x - m e) _ e(t )^)MM, 

da(x,t,e) dH^t—x + <f>(i),e) . d5(x — <j>(t), e) 
— -(Ti — h e(i) - 



dx g?£ d£ 

where denotes the derivative of u(., e) with respect to the hrst component. 

Substituting these relations into the left-hand side of the system (jl.ll) , we have the 
following 

du(x,t,e) du(x,t 7 e) da(x,t,e) 
___ +u{x>t>e) __ __ 

- uMt) d J^±mA +mR(x _ m*)-p(t)kt) W*-wu) 

dH u (-x + 0(t),e) 2 ^ dH u (-x + 0{t),e) 

- u ui — u-lH^-x + <p(t),e) — 



dH a (-x + 4>(t),e) d8(x-cf>(t),e) 
ai dC e(t) dl 



(4.2) 



and 



da(x,t,e) da(x,t,e) 2 du(x,t,e) 
___ + „(,., t , e )__ k ____ 



rfff g (-:c + ^(f),e) riff ff (-a: + 0(f),e) 
"001 37 Uio- x H u {-x + <p(t),e) — 

^- ^.e) + (t)ff + dS(x-J(tU) 
dt, a£ 

a lP (t)R(x - m ,e) d -^±mA +pmt)R{x - m e) W X -*®> e) 
k 2 u ^ dH u (-x + (f>(t),e) _ k% dR{x - <£(*), e) 



d£ r " ' d£ 

(4.3) 
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Now using the weak asymptotics relations from Lemma 3.1 in the relation (|4.2j) . 
we obtain 

du(x,t,e) , .du(x,t,e) da(x,t,e) 
___ +u(a:it>c )__ __ 

= {uKP(t) - uqUx - -u\ + ai}5(x - 4>{t)) + {^p 2 {t)uj - e(t)}5{x - 4>{t)) 
+ o v ,(l). 

Setting the coefficients of S and <5 in the above relation to be zero, we obtain 

Ux4>{t) - UqUx - ^llj + <Tl = 0, 

^p 2 (t)uj -e(t) = 0. 
The first equation above, when rewritten, gives 

*«> = %^ 

Substituting cb(t) from above in the relation (|4.3[) and observing that 

uM t)H u (- x + m t) ^MA + Uoe{t) d5 ( x e ) 



- e{t)<p(t) — = o V '(e), 

(the choice of c as in Section 3 helps us in getting this asymptotics) we have 

da(x,t,e) dcr(x,t,e) 7 du(x,t,e) 

-—— +<xM —^—-k ____ 

= {e(t) + crx4>{t) - 2 U i (T i - u o a i + k 2 u 1 }8(x - 4>(t)) + o v ,(l),e ->• 0. 
Setting the coefficient of S in the above relation to be zero, we obtain 
e(i) + <Ji4>(t) — 2 Ul(Jl ~ u ° ai k 2 ui = 0; 



which on simplification gives 



2 

e(t) = — - k 2 Ul . 



Ml 



Therefore we find that the smooth ansatz (|3.3p is a weak asymptotic solution 
provided pit), <fi(t), e(t) can be solved from the following equations 

4>{t) = 



r 2 



Ul 

^p 2 {t)uj -e(t)=0. 
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The above ordinary differential equations can be solved with the initial conditions 
0(0) = and e(0) = e° and we have 

<P(t) = [<7 \ , e (t)^(^l-k 2 Ul )t + e . (4.5) 

[u\ Ul 

Next substituting e(t) in the last equation of (|4.4[) . we can solve for p(t) taken in 
the form p(t) = pi(t) + ip2(t) and hence we have a weak asymptotic solution of 
the system (fTTjl . □ 

Since a 5 — shock wave type solution is the distributional limit of a weak 
asymptotic solution, we have from the previous theorem 

Theorem 4.2. For t E [0, co), the Cauchy problem (|l.l[) , (|3.2j) has a generalised 
8— shock wave type solution (|3.1I) with (f>(t) and e(t) given by the relations 

0(<) = [ ^ ] ~ [a] t, e(t) = - k 2 Ul )t + e°. (4.6) 

Remark 4.3. Let e° = 0. Now e(t) = in (|4.1|) would imply the existence of shock- 
wave solution {in the class of bounded variation functions) of the system (jl.ip . But 
that implies o\ = k 2 u\. Imposing Lax's admissibility condition, which in this case 
is u\ > 0, we get the shock curves 

Si : [a] = k[u], 

5*2 : [<t] = —k[u\. 

Also in this case, we have u^f- = lim e _;.o n{x, t, e) ^^^ ^ = ~ a i( u o + which 
again is the Volpert's product (the negative sign arises because of the convention on 
[a] ). Thus we recover the results proved in [B] for the shock-wave case (see (|1.3|l ). 

Remark 4.4. (Overcompressivity condition for 5— shock wave solutions) We recall 
that the overcompressivity condition (see PQ,[7]) for the S— shock wave solutions for 
a n x n system is 

Xk(vR) < <j>(t) < Xk(vL), k = l,..,n. 
Therefore for the system it takes the form 

[£] ~ M 
u - k < - < u + ui — k, 

[u] 

[£] - M 

uo + k < ^ < uo + u% + k. 

[u] 

The above relations on simplification yield the conditions: 

ui > 2k > 

and 

-(f-*X^<f-*. 

2 Ui 2 



Next we prove the existence of a weafc asymptotic solution for the system (|1.2p . 
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Theorem 4.5. For t £ [0, oo), the Cauchy problem (jl .2|) . ()3.2|) has a weak asymp- 
totic solution (13. 3[) wit/i <f>(t),e(t) and p(t) given by the relations 

N «i (4.7) 

^p 2 (t)io ~e(t) = 0, 

where ujq is a positive constant (defined in Section 3 ). 

Proof. Proceeding as in the proof of the Theorem 4.1, we find that the smooth 
ansatz (|3.3[) is a weafc asymptotic solution provided p(t),<f>(t),e(t) can be solved 
from the following equations 



(4.8) 



Ml 



ip 2 (t)w -e(t) =0. 



The ordinary differential equations for <fi(t) and e(t) can be solved with the initial 
conditions (f)(0) = and e(0) = e° and we have 

#t) = MzM t> e(<) = £f i + e o_ 

H mi 

Next substituting e(i) in the last equation of (|4.8[) . we can solve for taken in 
the form p(t) = p\(t) + ip%(f) an d hence we have a weak asymptotic solution of 
the system ([i~2|l . □ 



Therefore from the previous theorem, we have 

Theorem 4.6. For t G [0, oo), i/ie Cauchy problem (|1.2j) . (|3.2j) ftas a generalised 
5— shock wave type solution (|3.1I) «ra£/i </>(£) and e(i) given by the relations 

4>(t) = [ ^\~ [a] t, e(t) ^^-t + e°. (4.9) 
[u\ Ul 

Remark 4.7. T/ie overcompressivity assumption for the system (II. 2[) yields U\ > 
and-H± < ZL < «i 

2 ill 2 

Remark 4.8. I/ui > (follows from the overcompressivity condition above), then 
from equation (|4.7[) it follows that e(t) > 0. // m addition we have e° = 0, i/ien 

CT 2 

e(t) = -^f which is greater than zero for all t. In this case, it is sufficient to 
consider p(t) as a real-valued function only. 

Remark 4.9. From Remark 4-8 it also follows that we might have a singular con- 
centration developing in the solution of the system (|1.2[) even if we start with Rie- 
mann type initial data. 
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Remark 4.10. If we take <7o = o~i = in (13.2|l . then proceeding as in the proof of 
the Theorem 4-5 we obtain a generalised 5— shock wave type solution for the system 
(HI) of the form: 

u(x, t) = u + u\H(—x + <j)(i)), 

(4.10) 

a(x,t) =e°5(x-<f>{t)) y ' 

rail 

where <f)(t) = 

Remark 4.11. (Dependence of the solutions on k) 

From the structure of the generalised solutions for the systems , (|1.2[) obtained 
from Theorem 4.2 and Theorem 4-6, it's quite evident that as k tends to 0, the gen- 
eralised solution obtained for the system actually converges (in distributional 
limit) to that obtained for the system (|1.2[) . 

This observation therefore justifies our motivation to study the system (|1.2[) based 
upon the solutions of the system (|l.ip (letting k — > ). 
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